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Abstract 

Prom the recent measurement of parity nonconservation (PNC) in the Cs 
atom we have extracted the constant of the nuclear spin dependent electron- 
nucleon PNC interaction, k = 0.442(63); the anapole moment constant, K a = 
0.364(62); the strength of the PNC proton-nucleus potential, g p = 7.3 ± 
1.2 (exp.)±1.5 (theor.); the 7r-meson-nucleon interaction constant, f^ = h\ = 
[9.5±2.1 (exp.)±3.5 (theor.)] x 10 -7 ; and the strength of the neutron-nucleus 
potential, g n = — 1.7 ± 0.8 (exp.) ± 1.3 (theor.). 

PACS numbers: 11.30.Er, 21.10.Ky, 12.15.-y, 32.80.Ys 
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In the work |J the parity nonconserving (PNC) transition amplitude between the 6s and 
7s states of the 133 Cs atom has been precisely measured: 

E = -Im(£l pNC )//? = 1.5935(56)mV/cm. (1) 

They also observed the nuclear spin- dependent contribution 

Im(£l a )//3 = 0.077(ll)mV/cm. (2) 

This is a manifestation of parity violation in atomic nuclei and provides the first measure- 
ment of a nuclear anapole moment — an electromagnetic multipole violating the funda- 
mental symmetries of parity and charge conjugation invariance. The anapole moment was 
introduced by Zel'dovich pj just after the discovery of parity violation. He pointed out 
that a particle should have a parity- violating electromagnetic form factor, in addition to the 
usual electric and magnetic form factors. The first realistic example, the nuclear anapole 
moment, was considered in Ref. and calculated in Ref. [[|. In these works it was also 
demonstrated that atomic and molecular experiments could detect anapole moments. Sub- 
sequently, a number of experiments were performed in Paris, Boulder, Oxford, and Seattle 
H and some limits on the magnitude of the anapole moment were established. However, the 
first unambiguous detection of the nuclear anapole moment (14% accuracy) has just been 
completed |IJ. 

The existence of the anapole moment is due to parity nonconserving nuclear forces which 
create spin and magnetic moment helical structures inside the nucleus. (A detailed discussion 
of the spin helix produced by the weak interaction is contained in Ref. f6j). The wave function 
of the unpaired nucleon can be presented as 

i> = e ieoVr Vo, (3) 

i.e., the spin s = |cr is rotated around the vector r. Here the angle of rotation, 9r is 
proportional to the strength of the weak interaction [9 = ~^9P, see Eq. (|17D] and ipo is the 
unperturbed wave function. The correction to the electromagnetic currents due to this spin 
rotation has a toroidal structure. The toroidal electromagnetic current density j produces 
a magnetic field inside the torus like that inside a classical toroidal coil. In the limit of a 
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point-like nucleus the vector potential corresponding to this magnetic field can be presented 
as [g,| 

A = aS(r) 

a = — 7r f \(r)r 2 d 3 r = --^=—^——K a , (4) 
J JV ; ey/2 1(1 + 1) V ; 

where a is an anapole moment vector directed along the nuclear spin I, K — (I + 
|)(— l) I+1 l 2 ~ l {I is the orbital angular momentum of the external nucleon), and e is the 
electric charge of the proton. We separated the Fermi constant of the weak interaction 
(G) and introduced the dimensionless constant K a . The operator of the anapole moment, a 
(a = (^)|a|^)) is given by the following formula 0: 



7re 
m 



/i(r x <r) — ^(P r2 + r2 p) 



(5) 



where m is the mass of a nucleon, r and p are the position and momentum operators of the 
nucleon, fi is the nucleon magnetic moment in nuclear magnetons and q — (1) for a neutron 
(proton). The dominant contribution to the nuclear anapole is given by the spin current 
[the first term in Eq. @]. The contribution of the second term (the convection or orbital 
current contribution) is very small. Moreover, to a large extent it is canceled out by the 
contribution of the contact current (see Refs. []3],f|,|8|]). The only other sizeable contribution 
is due to the spin-orbit current considered in Ref. || and is about —20% of the dominant 
spin contribution. 

The interaction between atomic electrons and the magnetic field of the nuclear anapole 
produces a nuclear spin dependent PNC effect in atoms, which was first calculated in Ref. || 
and has been measured in Ref. (J). The PNC amplitudes for different hyperfine transitions 
were found to be different. This difference is produced by the magnetic interaction of the 
atomic electron and the anapole vector-potential A: 

V a = ea ■ A = ect ■ a.d(r) = ^j——-K a d(r). (6) 

Note that there are other mechanisms that produce (small) atomic effects similar to the 
anapole moment. This means that the atomic electron's interaction with the nucleus should 
actually be described by Eq. (||) with K a replaced by a new constant, k (more on this below). 
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Accurate atomic calculations of the PNC effect produced by the interaction have 
been done in Refs. |jTH|^T3[. The result of the many-body calculation in Ref. [IT] is very close 



to the semiempirical calculation in Ref. flQ |. The result of the Hartree-Fock calculation 
|fl"3|j differs by about 10% since it does not include many-body corrections. To reduce the 
theoretical error we calculate the ratio of the nuclear spin-dependent PNC amplitude to the 
main spin independent PNC amplitude. Using the most complete many-body calculation of 
the nuclear spin dependent PNC amplitude JTTJ , the calculation of the main PNC amplitude 
]TJj (which was done using the same method and computer codes) and the experimental 
data for different hyperfine transitions from Ref. |J we obtain the following equations: 

E(l + 0.05814k) = 1.6349(80) mV/cm (for the 6S F=4 -> 7S F=3 transition) 

E(l - 0.05148k) = 1.5576(77) mV/cm (for the 6S F=3 -> 7S F=i transition) (7) 



The solution to these equations is 



E = 1.5939(56)mV/cm 



(8) 



re = 0.442(63). 



(9) 



The calculated ratio of the nuclear spin dependent PNC amplitudes in ([?[) to the main PNC 
amplitude (E) is known very accurately, i.e., there is practically no theoretical error in the 
extracted value of re. This value of re contains three contributions: 



re — re a 



A" -1/2 1 + 1 



K 



re 2 + 



K 
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where K = 4 and / 
1.25(2 sin 2 6 W - \) « 



= | for 133 Cs, re a is the anapole moment contribution ([J), re 2 = 
-0.05 is the contribution of that part of the weak electron-nucleus 



interaction that depends on the nuclear spin (see, e.g., |6|||), and kq is the contribution of 
the combined action of the nuclear spin independent electron-nucleus weak interaction and 



the hyperfine interaction |15| (see also Refs. |Tq.p^ 

2.5 x lQ- 4 A 2/3 fi N = 0.017. 



1 a^ N 



(11) 
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Here Qw is the weak charge of the nucleus, a = e 2 = 1/137, = r v4 1//3 is the nuclear 
radius, and fi^ is the magnetic moment of the nucleus in nuclear magnetons (for 133 Cs 



/ijv = 2.58). The value of kq obtained in the more complete calculation in Ref. [|T7| is about 
1.5 times larger (as it contains some average radius of the nucleon distribution, R instead 
of Rn, as in the above equation), i.e., kq ~ 0.025. From the above results it follows that 
n a = 0.370(63). 

The Hamiltonian of the electron-nucleon interaction (|6p is presented for a point-like nu- 
cleus. However, a real nucleus has a finite size. Therefore, the "anapole moment" measured 
in the experiment |I] is in fact different than the anapole moment defined in (^). The 
"anapole moment" that was measured in the experiment can be defined as 



a = -7T J j(r)r 2 (l - Z 2 a 2 u{r)) d 3 r « (1 - 0.3ZV)a, (12) 

where u(r) ~ \ (^) (^) • For Cs Z 2 a 2 = 0.16. The interaction due to this "anapole 
moment" is just Eq. (|6]) with a replaced by a, i.e., the "anapole moment" is placed at 
the center of the nucleus. However, in the previous atomic calculations [|TUfTJ~5|] a different 
"regularization" prescription was used: 5(r) was replaced by a finite range function, 5(r) that 
has the shape of the nuclear density. The electron part of the anapole moment interaction 
([]) mainly mixes s\/2 and pi/2 electron orbitals. Using the electron wave functions inside 
the nucleus presented in Ref. || we have 

i . I ~ c*/ w i \ (V^lea • a<5(r)|VO 1 - 0.3Z 2 a 2 . , . ~. . . , . 

(^|ea ■ atf(r)|V P > = J _ 0A z* J = l-0.4Z^ ^ a|eQt ' a5(r)l ^" (13) 

This means that to accurately take into account the finite nuclear size the results of the 
atomic calculations of the anapole moment contribution pH|^r5| should be multiplied by 
(1 - 0.3Z 2 a 2 )/(l - 0AZ 2 a 2 ) w 1 + 0.1Z 2 a 2 = 1.016. Therefore the true value of « will be 
1.6% smaller than 0.37: 

n a = 0.364(62). (14) 



The value 0.36 has also been obtained in 19 



In Ref. 0] analytical and numerical calculations of n a have been done. The approximate 
analytical formula was obtained by using the wave function (|3p to calculate the mean value 
of the anapole moment operator (|]). The result is 

«« = ^—A 2/3 g P = 0.0H, (15) 
10 mr 



where fi is the magnetic moment of the external nucleon in nuclear magnetons and ro = 
1.2 fm. The more accurate numerical calculations ||,[S| in a Saxon- Woods potential with a 
spin-orbit correction give the following for 133 Cs: 



K a = 0.06g p . 



(16) 



Here g is the dimensionless strength constant in the weak nucleon-nucleus potential: 

G 

w = ■^7^[ (T ■ pp( r ) + pW 0- • p]> 



(17) 



where p(r) is the number density of core nucleons (g = g p for a proton). 

The proton-nucleus and neutron-nucleus constants can be expressed in terms of the 
meson-nucleon parity nonconserving interaction constants |)3],p0| (we use the notation of 
Ref. TO: 
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2.0 x 10 5 W n 



176^f / w - 19.5/1° - 4.7/iJ + 1.3^ - 11.3(/i° + /£) - 1.7// 



# n = 2.0 x 10 5 W/p 



-118^/, - 18.9/1° + 8Ahl - 1.3h 2 p - 12.8(/£ - /£) + Llfcj' 



(18) 



The parameters W p and are present in the above equation to take into account the 
nucleon-nucleon repulsion at small distances and the finite range of the true interaction 
potential. As in Ref. ||, we use the calculations of PNC for neutron and proton scattering 
on 4 He @, and take W p = 0.4 and W v = 0.16. Using the "best" values of the / and h 
constants listed in Ref. [^TJ (from here on we will refer to these as the DDH "best" values) 
one obtains g p = 4.5, g n = 0.2, and n a = 0.27. Note that this is a single particle shell-model 
value of the anapole moment constant. Shell-model calculations usually have an accuracy 
of about 30%. Thus, the agreement between the experimental value (0.364 ± 0.062) and the 
theoretical value (0.27) is as good as could be expected. (Moreover, it was shown in Ref. 
23| that the RPA corrections to the weak potential increase g p by 30%, thus K a could be 



increased to very close to the central experimental number of 0.364.) 

Comparing the measured value of K a ([14]) with the theoretical expression (|I6"D gives 



g p = 6± 1 (exp.). 



(19) 



We do not present here the theoretical error from the nuclear calculation of n a (about 30%). 



Now we can use the expression for g p in terms of the meson-nucleon interaction constants 
to find /„.. It was stated in the recent review that experiments give values of the p and 



uj weak constants very close to the DDH "best" values (these constants can be found from, 
e.g., p-p and p-a PNC experiments). The contribution of p and uj to g p is g p (p, uj) = 2. The 
main controversy is about the value of f n = h*- Comparison between flT9] ) and ( |T8"D gives 

/ w = {g p - 2) x 1.8 x 10~ 7 = (7 ± 2 (exp.)) x KT 7 . (20) 

We stress once more that the theoretical error in the nuclear calculation of n a is ignored 
here. Then, using this value of and the DDH "best" values of h p and in Eq. (|1^) we 
obtain 

g n = -0.38 x 10 7 f n + 1.9 = -0.9 ± 0.7 (exp.). (21) 

There are other nuclear calculations of K a [^,|l|,[T7],^(j] . Ref. [^5| contains a detailed 
calculation of the 7r-meson contribution to the anapole moment and Refs. Hl2|Jl7| include 
some configuration mixing effects. The most complete calculation of the anapole moment 
has been done in Ref. |26fl : they included all single-particle contributions (spin, spin-orbit, 
convection, and contact currents) and many-body corrections in the RPA approximation 
(e.g., the induced PNC interaction and the above mentioned RPA renormalization of the 
weak potential, which were considered in Refs. [T2^,[27]]). For comparison, it is convenient to 
present the result of their calculation in a form that stresses the role of g p . 

n a = 0.05{g p + 0.16g n - 0.07g pp - 0Mg np ), (22) 

where g pp and g np are the constants of the proton-proton and neutron-proton weak inter- 
actions; these are related to g p and g n by the formulae g p = (Z/A)g pp + (N/A)g pn and 
g n = [Zj A)g np + (N/A)g nn (see Ref. The authors of Ref. estimated the theoretical 

error in (E22j) as smaller than 20%. For the DDH "best" values of the meson-nucleon weak 



interaction constants we have g n = 0.2, g pp = 1.5, and g np = —2.2 |28| and so we obtain 



K a = 0.05(^- 0.05). (23) 

Thus, to an accuracy of ~ 1% K a is still proportional to g p . Comparing this with the 
experimental value of K a in Eq. flliD we obtain 



g p = 7.3 ± 1.2(exp.) ± 1.5(theor.). 



(24) 



Once again we can use the value of g p to find a value of /„■. Comparing the expression for 
g p (|18D with its numerical value (|24]) we obtain 



f v = hl= [9.5 ± 2.1 (exp.) ± 3.5 (theor.)] x 1(T 7 . (25) 

We increased the theoretical error here from 2.7 to 3.5 to take into account the uncertainty 
in the relation between g p and f n ([18]). As before, we use this value of /„- and the DDH 
"best" values of h p and h w in Eq. (|1^) and we obtain 

g n = -1.7 ± 0.8 (exp.) ± 1.3 (theor.). (26) 

We have presented two sets of estimates of g p , f w , and g n to give an indication of the 
possible spread of the results due to theoretical uncertainty. These two sets of results agree 
with each other to within their errors. In the abstract we presented values based on the 
more complete many-body calculations. 

Now we will compare our estimates of f n , Eqs. ( PD|) and (|2li|), with other estimates in the 
literature. There is no contradiction between these values of f n and the QCD calculations, 
which give f^ = h\ = 5-6 x 1(T 7 [0,0. The DDH "best" value of f n is / w = 4.6 x 1(T 7 . 



However, there are also smaller estimates of /„. in the literature, going down to the value 
\f n \ < 1.3 x 10~ 7 derived from a 18 F PNC measurement (see, e.g., the review ||31||). 



Note that there could also be a more exotic interpretation of the results of the k mea- 
surement: K2 may not be described by the standard electroweak theory and so may have a 
larger magnitude, thus implying a smaller value of k q [see Eq. fllOD l, and hence f n . However, 
such an explanation would be very improbable since the results of measurements of atomic 
weak charges and PNC in deep inelastic electron-nucleon scattering agree with the standard 
model. To clear this question it would be interesting to measure the anapole moment of 
the 207 Pb nucleus, which contains an external neutron. The constant g n contains f n with a 
negative sign in this case [see Eq. (|IBD]. 

Just before the submission of this paper it was brought to our attention that an analysis 
of nucleon weak interactions, based on the experiment [0, has also been done in the recent 
work p2|. 



One of us (VVF) is grateful to LB. Khriplovich for useful comments. This work was 
supported by the Australian Research Council. 
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